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The flow equation approach investigated by Wegner et al. is applied to an unbounded 

Hamiltonian system with a generalization. We show that a well-known quantized complex energy 
eigenvalues which is related to decay widths can be given with this approach. 
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To investigate a quantum system, a diagonalization of 
a given Hamiltonian is one of mighty approaches. Eigen- 
values and eigenstates which can be obtained with the 
diagonalized representation are most important clues to 
understanding of the physical properties of the system. 
But unfortunately, however there are a few explicitly 
solved exceptions, an analytically diagonalization of a 
given Hamiltonian is hard problem and various schemes 
have been deviced. Among of them, a novel approach 
to this problem was proposed by Wegner jjj. A way to 
construct of a generator, which generates one parameter 
unitary group which drives a given initial Hamiltonian 
to diagonalized (or almost diagonalized) final one, was 
deviced. This approach is called flow equation approach 
and it has already been applied to various physical mod- 
els, especially condensed matter physics and high energy 
physics with some techniqal improvements 

On the other hand, there is interesting argument on 
the diagonalization, which is the diagonal representation 
of a Hamiltonian unbounded from below with complex 
eigenvalues. The imaginary part of the eigenvalues is re- 
lated to the life-time of the quasi-stationary state, and it 
plays a good role to investigate a resonance phenomena 
H |To| , however the eigenfunctions belonging to the com- 
plex eigenvalues can not be interpreted in a conventional 
sense. A mathematical treatment of the quasi-stationary 
approach has been also studied with a rigged Hilbcrt 
space (or Gel'fand triplet) formulation of quantum me- 
chanics [|]-[l0| . As pointed out by Bohm etal || , such for- 
mulation allows for the description of the intrinsic irre- 
versible processes in quantum mechanics, and these stud- 
ies would be important in order to understand the fun- 
damental relation between macroscopic irreversible phe- 
nomena and microscopic reversible theory. 

Our interest is an investigation or generalization 
of Wegner's flow equation to the unbounded system 
problems. In this letter, we consider one-dimensional 
quadratic interaction model (Q), which becomes un- 
bounded below in cases of |Ao/wo| > 1/2. We show that 
the flow equation approach can be generalized to such sit- 



uation where the original approach does not work, and 
this treatment provides a good scheme to study such deep 
features as a time arrow in quantum mechanics which 
were discussed in fls|-|iof . 

We shall consider a quadratic interaction system with 
the following Hamiltonian 



H 



ujq a 1 a + Ao (aJ + a 2 ^j + 



vo 



(1) 



where uj , X and Vg are given constant of real numbers. 

In order to explain the wegner's approach and our in- 
terest, first, let us apply the original flow equation ap- 
proach to this Hamiltonian (|l]). The flow equation is a 
continuous unitary transformation of the Hamiltonian, 
which is written in the following differential form; 



dH(l) 
all 



= [r,(l),H(l)], r l (l) = [H d (l),H(l)} (2a) 



H(l) = U\l)H(0)U(l), -U\l)=r,(l)UHl) (2b) 

with the initial condition H(Q) — H (original (given) 
Hamiltonian). Hd(l) is the diagonal portion of H (I). Ap- 
plying to the Hamiltonian (Q) , we obtain the following set 
of differential equations by coefficient matching: 



dcu(l) 

dl 

d\{l) 

dl 
dv(l) 

dl 



= -Ww(l) A(0 2 
= -4w(0 2 A(0 
= -8 A(0 2 , 



(3a) 
(3b) 
(3c) 



where 



H(l) = tj a\l)a(l) + Ao U\lf + a{l) 



){l)a)a + \{l) (a) 2 + a 2 ) +v(l) (4a) 



H d {l)=u(l)a^a + v(l), 



(4b) 
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and 



and 



t(/) = U\l)a)U{l), a(l) = U\l)aU(l) (4c) 

with cj(0) = wo, A(O) = Ao, v(Q) = Vq. One can easily 
find that the solution of these equations has satisfied 

lu{1) 2 - (2 A(/)) 2 =u>l-{2 A ) 2 = const. (5a) 
v(1)=v + -(uj(1)-uj ). (5b) 

Notice that there are two different type solutions whose 
behaviors are qualitatively different from each other 
(Fig.0). 



h= -1/2(0 




X= 1/2 co 



FIG. 1. Schematical illustration of the unitary flow. The 
filled circles and crosses represent the fixed points (n) and (fiy. 



They can be distinguished by their initial points, 
(i) |Ao/wo| < 1/2 case 



u(oc) -> w y 1-4 (J^- ) . 
A(oo) -> 0, 

1 



u(oo) -> w + ^ - w ) 

(ii) (Ao/wol > 1/2 (unbounded) case 
a>(oo) — > 0, 



A(oo) -*■ A y 1 - 
v(oo) -> w - 2 W o- 



1 / LUq 



4 VAo/ ' 



(6a) 
(6b) 
(6c) 

(7a) 
(7b) 
(7c) 



Notice that the diagonal representation (A(oo) = 0) is 
obtained in the case of (i) (|Ao/wo| < 1/2). This diagonal 
representation (|J) can be understood in terms of famous 
Bogoliubov transformation, i.e. 



H = atfb + K + vo, [6,6 f ] = l 



where 
Q = 



U 



sinh 2 9 + cosh 2 9 ' 



uiq sinh 9 
sinh 2 9 + cosh 2 ' 



(8a) 



(8b) 



a = 6 cosh — sinh#, = cosh — 6 sinh 9, (9a) 

with which satisfies 

Aq sinh 9 cosh 9 



sinh 2 + cosh 2 9 



(9b) 



One can easily check that relations (||) are completely 
equivalent to M), i.e., 



u)(oo) = f2, v(oo) — K + t>Q. 



(10) 



The eigenvalues of (Q) can be easily obtained from (Q), 
that is 



£"„ = Q n + /t + Vq (n is integer). 



(11) 



As it should be, the original flow equation approach 
works in this case (i). 

In not only this case but also general cases where we 
can obtain a diagonal representation of a given Hamilto- 
nian with this approach, the flow equation approach or 
the unitary transformation U(oo) can be interpretated as 
a transformation of the basis from one which diagonal- 
izes Hd(0) to another which diagonalize H(0). We can 
say this original approach works when both of the basis 
are in a same function space because U(oo) is unitary. 
On the other hand, in the case of (ii) it is clear that the 
eigenstates of the Hamiltonian is only scattering states 
and they can not be in a same function space with eigen- 
states of Hd(0) (boundary state). Therefore it is rather 
reasonable that the original approach does not work, as 
we saw in (Q). 

It would be obvious that some generalization to larger 
class flow should be required to obtain the diagonal rep- 
resentation of the Hamiltonian. The key should be to 
remove the unitarity from the flow equation approach 
and that it can be realized in the following procedure. 

Generally, flows with the differential equation of the 
first-order like equation (^) does not have any crossing 
points with each other, except for the attractor points. 
Thus we would expect that some different flows whose at- 
tractor points are common will be related to each other. 
A diagonal representation of an unbounded Hamiltonian 
is also obtained with some new flow which has same at- 
tractor points with the original (unitary) one, however 
this new flow would not represent an unitary transfor- 
mation, unlike the original flow. General solutions of 
Eq. (||) which has a attractor point (0) can be written as 

u>(l) = 2A(oo) (cos<pi sinh£; + ismipi cosh£/) , (12a) 
A(^) = A(oo) (cos ipi cosh£; + ismtpt sinh£/) , (12b) 



v(l) = v + - MO - wo) , 
where 



(12c) 
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dfi 



dl* 1 



-8A(oo) cos 2(fi sinh 2£z, 
-8A(oo) sin 2ip t cosh 2£/ 



with 



Coo — ► 0, (fi 



0. 



(13a) 
(13b) 

(13c) 



It should be noticed that the real (unitary) flow is one of 
the solutions represented with (p^), i.e., tpi = 0, and there 
are two other separable solutions i.e., ipi > ((fi — ► +0 
as / — > oo) and < (<^; — > — as / — » oo). which are 
originated from unstable fixed-points 

(w,A) = (±2*A(oo),0), (ar(£ )¥ >) = (0,±~)). (14) 




Re?i 



Re co (oo^o ) 

FIG. 2. Schematical illustration of the generalized flows 
# (0) (Z) and H (+) (Z). 

There are three kinds of flows which have a common 
attractor point (0, A(oo)) but different "initial (I — > 
— oo)" points, (2iA(oo),0), (— 2iA(oo), 0), and "some- 
where on the real flow plane" . Corresponding to them, 
we can introduce three Hamiltonian- flows, i.e., H^\l) 
and HW(l) as 

H {}) {I) = {l)a)a + \® (Z) (a) 2 + a 2 ) + u« (i) (15) 

(i = {±,o», 

with 

w (± )(-oo) -► ±2iA(oc), A (± )(-oc) -> 0, 

(v (±) (-oo) -> v{oo) ±zA(oo)) (16a) 

and 

J°\0) = loo, A(°)(0) = A , («(°)(0) = «o) . (16b) 

It is obvious that these Hamiltonian-flows are attracted 
to the common fixed point (0) as I — > oo, i.e., 

w (j) (oo) -> 0, A (j) (oo) -> A(oo), v u) {oo) -> u(oo). 

(16c) 



Precisely speaking, ( |l3| ) and (16) are not sufficient to de- 
fine H^\l) uniquely. They depend on direction of an 
"infinitesimal shift" from the fixed points. Giving this 
direction is giving some complex number V5oo/£oo- For 
example, <^oo/£oo = corresponds to H^°'(l). 

The Hamiltonian flow ff'°)(/) is the unitary flow de- 
scribed in (0), and that not H^(l) but H^(l) give us 
the diagonal representation of the original Hamiltonian 
([j]) with complex eigenvalues (at I — > — oo), i.e. 

ff (±) (-oo) = ±2iA(oo) a)a + v (±) (-oo) 




+ v(oo) 



(17) 



We can immediately obtain the quantized complex eigen- 
values 



- ) -!- r( x i 



(n is integer), (18) 



and this is consistent with well-known result (for exam- 
ple, see ^,0). The symmetry in ( |l8| ) is related to sym- 
metric phenomena in unstable system, that is, "decaying 
process for the positive time" and "growing process for 
the negative time", as was pointed out in ||[l0|]. On 
should notice that the infinitesimal shift from the attrac- 
tor point (0) decides the arrow of time of the phenom- 
ena. This means that the structure of the generalized 
flow should be directly connected to the deep features 
of an irreversible process which described with quantum 
mechanics. 

Let us mention to applications of our approach to other 
models. 

First, this approach seems possible to apply easily to 
other models, however we have to be careful of the fol- 
lowing fact. In general, when our approach is applied to 
some systems, a truncation schemes to get a closed set 
of differential equations would be often required, since 
the higher interactions will be generated with the "evo- 
lution" [/(/). In such cases, we can adopt some trunca- 
tion schemes which have already been proposed for the 
ordinary flow equation approach ]3|,|| in which we can 
consider non-perturbative corrections. 

Second, in general cases it is almost impossible to solve 
the set of equations analytically because they becomes 
highly nonlinear even after the trancation. Therefore, 
sometimes, we have to study generalized flow equations 
numerically. In the numerical studies, the following pro- 
cedure to find an unstable fixed point (which gives a di- 
agonal representation of unbounded Hamiltonian) would 
be useful. 

1. Find a stable fixed point (like urn) with the forward 
(I — > oo) evolution with the unitary flow from orig- 
inal (given) Hamiltonian. This process should be 
always easy because the point is a attractor point. 
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2. After finding the fixed point, compute the backward 
evolution from the "initial" point which is slightly 
(infmitesimally) shifted to some direction from the 
attractor point. 

The flow will be "attracted" to the unstable point as 
I — > — oo, and we should easily find the unstable fixed 
points as "attractors" . There must be such symmetric 
two unstable points which depend on the direction of 
the "initial" infinitesimal shift. (This corresponds to two 
Hamiltonian flows ijW(Z) for (Q). ) 

Finally, let us mention to an application of our ap- 
proach to a so-called open system problems. As well 
known, there are two kinds of discussion about the time 
arrow in quantum mechanics [||. One is called the in- 
trinsic irreversibility and the other is called the extrinsic 
irreversibility. The decaying process of state with the 
Hamiltonian (]]]) is thought as an example of the former, 
and the decay with an interaction with environment is 
thought as the latter or as a dynamics of the open system 
HI!. Dealing with the open system is another stan- 
dard approach to introduce the time arrow to quantum 
mechanics. It should be interesting problem to apply our 
approach to open systems. Actually spin-boson model 
which is a typical model to discuss the extrinsic decaying 
process Jl2| , was studied within an (improved but) uni- 
tary flow equation approach by Mielke etal Q , and they 
succeeded in computing the energy shift of spin system 
which is due to the interaction with the bosonic environ- 
ment. But decay constants was not obtained directly and 
the dissipative property of the system seems hard to be 
understood from their analysis. We think our approach 
can be used to reexamine this system and the dissipative 
properties of this system can be made clearest. More- 
over it will be helpful to understand the relation between 
two different treatment on the irreversibility in quantum 
mechanics. 
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